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permits the use of radiosonde data to provide model 
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A MODEL ATMOSPHERE FOR EARTH RESOURCES APPLICATIONS 

By David E. Pitts and Kirby D, Kyle 
Manned Spacecraft Center 

SUMMARY 


All earth resources remote- sensing techniques are affected by the atmosphere 
lying between the target and the sensor. The computer program presented in this re- 
port offers a method of numerical use of radiosonde data so that atmospheric effects 
may be assessed and possibly removed from the signal. 


INTRODUCTION 


The objectives of the NASA Earth Resources Program are to determine the per- 
formance capabilities of various sensors, to discover signature criteria of resources, 
and to develop new sensors and systems that will eventually enable management of 
earth resources. To accomplish these objectives, certain absolutes which may be used 
to evaluate sensing systems and techniques must be established. The laboratory usu- 
ally offers the best testing environment, but the type of target, the conditions of the 
path of the signal, and other testing parameters are limited. In general, the labora- 
tory is so restrictive that a successful laboratory test of a remote sensor is necessary 
but not sufficient to ensure proper operation of the sensor in an application. There- 
fore, much of the testing is performed in the same environment in which the instrument 
is expected to operate. Testing under such conditions requires that the data concern- 
ing the environment between the instrument platform (e. g. , an aircraft or a spacecraft) 
and the target be as accurate as possible. Thus, determination of the "ground truth" 
and description of the state of the atmosphere in the path of the electromagnetic signal 
are necessary. 

Remote- sensing techniques are affected by the atmosphere lying between the tar- 
get and the sensor. The amount of noise introduced into the signal by the interaction 
between the atmosphere and the signal depends upon the type of sensor, the wavelength 
employed, and the meteorological conditions prevailing at the time of the experiment. 
Since the NASA Earth Resources Program remote- sensing effort is in a developmental 
stage, the effects of this interaction are presently being determined, and hopefully, 
the model atmosphere for earth resources applications, presented in this paper, will 
facilitate analyses of such effects. 

The computer subprogram set presented in this paper offers a self-consistent 
method for numerically calculating the state of the atmosphere based on radiosonde 



data given in terms of significant levels of pressure, temperature, and temperature- 
dewpoint depression. After data from the radiosonde closest to an aircraft or space- 
craft remote -sensing target have been obtained and after these data have been inserted 
into the computer subprogram set, a programer has almost any desirable atmospheric 
parameter available for use in his computer programs. In particular, the subprogram 
set described in this paper makes available all the necessary quantities for calculation 
of infrared and microwave absorption or refraction, or both. However, no attempt has 
been made in this paper to include atmospheric absorption calculations in the model at- 
mosphere; only the basic atmospheric data necessary for the previously mentioned cal- 
culations are provided. 

The model atmosphere was written in the FORTRAN V computer language for the 
Univac 1108 computer. However, the program is also compatible with Control Data 
Corporation and IBM FORTRAN IV compilers. Copies of the computer cards are avail- 
able upon request from David E. Pitts, TF8, Manned Spacecraft Center, Houston, 

Texas 77058. 


SYMBOLS 


^ sound 

computer symbol 

d 

computer symbol 

e 

computer symbol 

e s 

computer symbol 

f w 

computer symbol 
vapor mixture, 

g 

computer symbol 

g 0 

surface gravity, 

H 

computer symbol 

H 

a 

computer symbol 

H 

b 

computer symbol 

H 

P 

computer symbol 

H 

P 

computer symbol 


ANS(4), speed of sound, m/sec 

DUM/D1, increment of the slant path from r to r\ cm 

ANS(19), E(X), water-vapor pressure, mbar 

ANS(20), E(X), saturation water-vapor pressure, mbar 

F(P, X), correction for the departure of the air and water- 
from ideal-gas law 

2 

ANS(5), acceleration caused by gravity, f(Z), cm/sec 
2 

g at R , cm/sec 
° e 

H(I), geopotential altitude, m 

HA, HLOW, geopotential altitude at A, where H a < H^, m. 

HB, geopotential altitude at B, m 
ANS(15), pressure scale height, km 
ANS(16), density scale height, km 
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M mass percentage of the ith constituent 

i 

m computer symbol ANS (7), molecular weight of the atmosphere, g/g-mole 

m, Q molecular weight at H^, g/g-mole 

computer symbol XMG, molecular weight of the dry atmosphere, g/g-mole 

m computer symbol XMO, molecular weight at the surface, g/g-mole 

o 

m molecular weight of water, g/g-mole 

w 


n’ 

computer 

n" 

computer 

n STP 

computer 

n(Z) 

computer 
and P 

P 

computer 

P 

a 

computer 

P v 

b 

computer 

q 

computer 

q 

H s 

computer 

R 

computer 

R 

e 

computer 

Rel 

computer 

Pl x 

computer 

Ry 

computer 

R z 

computer 


symbol XN2, refractive index at r' + (1/2)AZ 
symbol XN1, refractive index at r" + (1/2) AZ 
symbol ANS (17), refractive index of air at STP 

symbol ANS(18), refractive index of air as a function of X, T, 

symbol ANS(l), atmospheric pressure, mbar 
symbol PLOW, atmospheric pressure at H , mbar 

symbol PHIGH, atmospheric pressure at H^, mbar 

symbol ANS (13), specific humidity, g/kg 

symbol ANS (14), specific humidity at saturation, g/kg 
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symbol RO, universal gas constant, 8. 31432 x 10 ergs/(mole °K) 
symbol RE, mean radius of the earth, 6371. 299 km 


symbol ANS (12), relative humidity, percent 


symbol XS-XL, X- component of (r ' - r^J, 




symbol YS-YL, Y- component of ^r ' 


T 

sp 1 


km 


km 


symbol HS-HL, Z-component of ' - r^ j, 


km 
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r 


computer symbol ANS(IO), mixing ratio of the water in the atmosphere, g/kg 
r' computer symbol S2, distance to shell Z + AZ on the refracted path, km 

r" computer symbol SI, distance to shell Z on the refracted path, km 

r^ distance from the center of the earth to a target, km 

r g computer symbol ANS(ll), mixing ratio required for the saturation of water 

in the atmosphere, g/kg 

r ' distance from the center of the earth to a spacecraft, km 

S computer symbol S, Sutherland's constant, 110.4 °K 

s distance 

T computer symbol ANS(2), kinetic atmospheric temperature, °K 

T* computer symbol ANS(6), virtual temperature, °K 

T computer symbol T( ), temperature at H , °K 

R cl 4 

T * computer symbol TVLOW, virtual temperature at H , °K 
computer symbol ANS(2), temperature at H^, °K 
T^* computer symbol TVHIGH, virtual temperature at H^, °K 

dewpoint temperature, °K 

T , computer symbol TD( ), dewpoint temperature at H , °K 

b computer symbol ANS(9), dewpoint temperature at H^, °K 
TV molecular scale temperature, °K 

TOS computer symbol TOS, angle betv/een r' and r ', rad 

1 sp 

t time 

VV identifier of the significant- level data set of radiosonde code 

Z computer symbol Z, geometric altitude, km 
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computer symbol ZL, altitude of a target above the earth, km 

Z computer symbol ZS, altitude of a spacecraft above the earth, km 

/3 computer symbol BETA, 1. 458 x 10"® kg/ (sec °K m) 

y ratio of specific heats 

£ computer symbol PHI, angle from the zenith down to the tangent to the path 

at the target, rad 

computer symbol C(3), distance upward from a local station to a spacecraft, 
rad 

tj" computer symbol C( 2), distance eastward from a local station to a space- 

craft, rad 

0 j computer symbol THETAL, target longitude, input card, deg (internally, rad) 

0 computer symbol THETAS, spacecraft longitude, input card, deg (internally, 

sp rad) 

X computer symbol XLAMDA, wavelength, microns 

H computer symbol ANS(8), coefficient of viscosity, kg/ (msec) 

£ computer symbol SUM1, dummy variable, rad 

ij" computer symbol C(l), distance southward from a local station to a space- 

craft, rad 

3 

p computer symbol ANS(3), atmospheric density, g/cm 

3 

p^ density of dry air, g/cm 

3 

p density of water vapor, g/cm 

4>' computer symbol PHIPR, angle between r' and the path of the ray after 

refraction, rad 

<p" computer symbol PHI, angle between r" and d, rad 

computer symbol PHIL, target latitude, input card, deg (internally, rad) 
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computer symbol PHIS, spacecraft latitude, input card, deg (internally, rad) 


xfy computer symbol PSI, angle between r' and d, rad 

MODEL ATMOSPHERES 


Model atmospheres for earth resources applications may be described as one of 
three types: preflight, flight, and postflight. Preflight model atmospheres include 
those which have been developed from aerospace flight- support models (refs. 1 and 2) 
and statistical models of cloud cover over the earth (ref. 3). The last of these 
indicates the probability of success on spacecraft- or aircraft-borne photographic 
tnissions for earth resources applications. 

Flight model atmospheres are calculated from sounding-type remote- sensing de- 
vices aboard spacecraft or aircraft. Flight model atmospheres are not presently well 
developed, but when they are well developed, they will represent the ultimate in knowl- 
edge of the "air truth” until special-purpose instruments that will perform atmospheric 
noise extraction in real time are developed. 

Postflight model atmospheres are based upon standard meteorological soundings 
and are used to assist in the development of flight model atmospheres. These post- 
flight model atmospheres may be described as predictive and nonpredictive. 

Predictive postflight model atmospheres use equations of motion, thermodynam- 
ics, and continuity and standard meteorological soundings to predict (in time and space) 
the state of the atmosphere near the target for a remote sensor mounted on an instru- 
ment platform. This type of model atmosphere is not presently well developed. Non- 
predictive postflight model atmospheres offer a self-consistent method of calculating a 
model atmosphere at the position of a radiosonde which may be located near the experi- 
ment platform. The subprogram model atmosphere set discussed in this paper has the 
capability of performing either as a nonpredictive postflight model atmosphere or as a 
preflight model atmosphere, depending on the form of the input data. 


EQUATIONS FOR THE MODEL ATMOSPHERE 


The model atmosphere may generally be considered to be in a state of quasi- 
static equilibrium. That is, when the equations of motion, thermodynamics, and con- 
tinuity are scaled and when closed sets are found, the large-scale (i. e. , the first 
order) vertical- component solution will show that, except near clouds with high-velocity 
updrafts, the hydrostatic equation 



(1) 
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applies well. In equation (1), P is atmospheric pressure, Z is geometric altitude, 
p is atmospheric density, and g is the acceleration caused by gravity. At pressures 
and temperatures experienced in the atmosphere of the earth, the ideal-gas law is usu- 
ally accurate to within 1 percent. The equation of state 


- ^ >m 
p “ RT 


( 2 ) 


is a form of the ideal-gas law, where m is the molecular weight of the atmosphere, 

R is the universal gas constant, and T is the kinetic atmospheric temperature. 

With certain reasonable and valid assumptions, the proper combination of the 
hydrostatic equation (eq. (1)) and the ideal-gas law (eq. (2)) results in equations (3) and 
(4), which are derived in detail in reference 4. If 3T*/ 3H £ 0, where T* is virtual 
temperature and H is the geopotential altitude, then 





g o m d /[R(3TV3H)] 


( 3 ) 


and if 3T*/ 3H = 0, then 


P, = P exp 
b a 1 


S„ m d( H b - H a)l 


RT * 

cl 


(4) 


In equations (3) and (4), P, is the atmospheric pressure at H, , P is the at- 

u d a 

mospheric pressure at H a , T^* is the virtual temperature at H^, T' a * is the virtual 

temperature at H , g is the surface gravity, m, is the molecular weight of the dry 

atmosphere, H is the geopotential altitude at A, and H, is the geopotential altitude 

at B. In the upper atmosphere, a fictitious temperature designated as molecular 
scale temperature T m is defined in order to include variations in molecular weight 

(caused by molecular dissociation) and temperature in one variable. 

m 

rp _ rp 0 

m m 


( 5 ) 



where m Q is the molecular weight at the surface. Similarly, in the lower atmosphere, 

a quantity designated as virtual temperature T* is defined in order to include varia- 
tions in molecular weight (caused by water vapor) and temperature in one variable. 


m 


T* = T 


m 


(6) 


Therefore, T* and T may be used interchangeably in equations (3) and (4); this 

fact enables the use of equations (3) and (4), which were derived for planetary atmos- 
pheres in reference 4. 

As shown in appendix A, the proper combination of the equation of the state of 
dry air, the equation of the state of moist air, and equation (6) gives the exact expres- 
sion of T* as a function of temperature, pressure, and water-vapor pressure. 


T* = 


1 - 0. 37803 


f e 
w 


( 7 ) 


where f w is the correction factor for the departure of the air and water-vapor mix- 
ture (from the ideal-gas law) and e is water-vapor pressure. Equations (3) and (4), 
which are the fundamental equations of subroutine MODATM calculations, are used in 
different forms to find the altitude of the significant levels and to find the pressure at 
a level between significant levels. 


Subroutine MODATM 

When atmospheric data at a particular altitude are desired, either geometric 
altitude is used as the calling variable, or pressure is used as the calling variable and 
a corresponding geometric altitude is calculated by using equations (3) and (4). Geo- 
potential altitude H is calculated by 


H = 


Z ( R e) 

R + Z 
e 


( 8 ) 


where R g is the mean radius of the earth. Geopotential altitude is then used to cal- 
culate temperature, virtual temperature, and molecular weight. 
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Temperature is calculated by 


T. = T 
b a 


+ S(H„ 


3H b 


H 


(9) 


where is the temperature at and T a is the temperature at H . Virtual tem 
perature is calculated by 


T 1 * = T * + 
b a 


3T* 

3H 


( H b - H a) 


( 10 ) 


Molecular weight is calculated by 


m ,T 

™ _ d b 

m b - T * 


( 11 ) 


where m^ is the molecular weight at H^. 

When P and T* are known, a form of the equation of state (eq. (2)) 

Pm, 

P = "RT* 


is used to calculate density. Then, additional quantities related to altitude, pressure, 
density, molecular weight, temperature, and virtual temperature are calculated. The 
equations for the speed of sound C goun( j, acceleration of gravity g, coefficient of vis- 
cosity p, saturation mixing ratio r , saturation specific humidity q , pressure scale 

s s 

height Hp, and density scale height are as follows: 


C = / RT* 
sound m^ 


g g o ( r b + z) 
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(15) 


H = 


/3T 3 / 2 

T + S 


0. 62197f e 
w s 

/P - f ej 
\ w s) 


(16) 


0. 62197f e 
w s 

q s ~rp - 0 . 37803f e X 
\ w s) 


(17) 


H - RT * 
P" m d g 


(18) 


H p " T 1 7 3T*\ 

H + T* v az / 
p 


(19) 


- fi 

where y is the ratio of specific heats, /3 is 1.458 x 10 , S is Sutherland's constant, 

and e is the saturation water-vapor pressure. Equations (13), (15), (18), and (19) 
s 

are derived in reference 1, equation (14) is derived in reference 4, and equations (16) 
and (17) are derived in reference 5. The f w -factor is calculated by a function subpro- 
gram simulating tables 89 and 90 given in reference 6. 

For calculations of variables describing the amount of water vapor in the atmos- 
phere, dewpoint temperature is calculated as follows: 


3T ■ 


T , , — T . + 

d, b d, a 


-==(H. - H ) 
3H \ b a/ 


( 20 ) 


where b is the dewpoint temperature at H^, and a is the dewpoint tempera- 
ture at H . The equilibrium vapor pressure over a plane surface of water (ref. 6) is 
then calculated. 
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. e = 1013. 246 X 10 


-7.90298^-1. 0+(373. l^TjJ+5. 02808 log 10 (373. l^T^)-!. 3818x10' 


10 


11.344[l.0-(T d /373. 16) 


1.0j+3. 1328xl0' 3 ho 


-3. 4914 [-1. 0+ .'373. 16/T d )j 


(21) 


The formula for the vapor pressure over ice (ref. 6) may also be used. 


-9.09718 

e = 6. 1071 x 10 


£ 


1. 0+(273. 16/T d )]-3. 56654 log 10 (273. 16/Tj)+0. 876793]!. 0-(Tj/2?3. 16)] 


(22) 


The choice of the temperature ranges during which each of the previously mentioned equations for e is used is 
determined by the programer (function E(X)). As presently set up, only equation (21) is used. Equations (21) and 
(22) are used for calculating e by using T in place of T . 

S Cf 

With the previously discussed basic quantities available, the remaining atmospheric quantities may be cal- 
culated. The equations for the mixing ratio r, relative humidity Eel, specific humidity q, refractive index 
ng^p (in wavelength), and refractive index n(Z) (in P, T, and wavelength) are as follows (ref. 5) : 


0. 62197f e 
w 

= ( p ' f w e ) 


(23) 


Rel = — x 100 
r 

s 


(24) 


0. 62197f e 
w 

q " JP - 0. 37803f^eJ 


(25) 



Fox’ the infrared region (ref. 7) 


^TP 


= 1 + 10“ 8 /6432=8 + 


2949810.0 25540 


146 


1 


41 - 


and 


A 288. 15\ 

1 + 273. 16 \ P 


n(Z) = l +( n STp -l)— ^ 

V + 273. 16/ 


1013.25 


(26) 


(27) 


where A is wavelength. If the wavelength is in the microwave region (A > 12 500 mi- 
crons, i. e. , A > 1. 25 centimeters), then 


n(Z) = 1.0 + 


1.0 x 10' 6 



373000. 0 



(28) 


as shown in reference 8. 

The input variables of MODATM are included in the calling argument, and all 
output variables (i. e. , the variables calculated by equations (3) to (28)) are stored in 
a "common block" in the array ANS. Detailed instructions on the use of subroutine 
MODATM are included in comment cards. For data-card information, see the discus- 
sion on subroutine INPUT in this report. 


Subroutine INPUT 

The purpose of subroutine INPUT is to read the input data cards necessary to set 
up the significant levels of various atmospheric parameters (i. e. , altitude, pressure, 
temperature, and dewpoint temperature) for subroutine MODATM. Subroutine INPUT 
is initiated by MODATM whenever pressure (i. e. , ANS(l)) is set equal to a number 
which is less than zero, and because of this fact, many sets of radiosonde data may be 
used successively, but not concurrently. 

The input data may be of the form given in the significant levels (i. e. , VV) of 
pressure, temperature, and temperature-dewpoint depression for a radiosonde. 
Table I shows an example of radiosonde data and the key to the radiosonde code. 
Table II gives the input data cards for the example shown in table I. 

Subroutine INPUT is also constructed to accept input data other than radiosonde 
code VV. If the first data card encountered is blank, then each of the next data cards 
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will be read in uneoded form (i. e. , as altitude, temperature, and relative humidity). 
An example of the input data cards necessary to set up the 15° N annual model (ref. 2) 
is included in table III. 

Levels of possible condensation are indicated by the word "condensation” 
in the print-out of the significant levels. This occurrence is determined by 

T - T f| < 2° K at 1500 meters and T - < 8° K at 9000 meters, which is expressed 

by the approximate expression 


T - T d < 1. 0 + 0. 000777H (meters) 


(29) 


Subroutine REFRAC 

Subroutine REFRAC is included to assist in making refracted path calculations 
throughout the atmosphere. The basic equations are developed (ref. 9) from Snell’s 
law 


n’ sin 0 ’ = n" sin 0 (30) 


and from the law of sines 


sin 0" _ sin 0 
r' - r" 


(31) 


as shown in figure 1. In equations (30) and (31), n' is the refractive index at 
r' + (1/2) AZ, 0' is the angle between r' and the path of the ray after refraction, 
n" is the refractive index at r" + (1/2)AZ, 0 is the angle between r' and d, 0" 
is the angle between r’ ' and d, r’ is the distance to shell Z + AZ on the refracted 
path, and r" is the distance to shell Z on the refracted path. 

The combination of equations (30) and (31) gives 


. -l/n”r” sin 0"\ 
* sem { Hv ) 


(32) 


and 


0 = sin 


- 1 (^_ 


sin 0" 
r’ 


(33) 



Thus, by using known values for r", r', A, and 0" and by initiating MODATM to 
obtain values for n" and n\ the angles 0' and \f/ are calculated. If a continuous 
path is desired, 0” should be set equal to 0', and r" and r' should be incre- 
mented. Then, subroutine REFRAC should be called again. 

Slant-path calculations are also made available by using the law of sines to cal- 
culate the increment d of the slant path from r to r' as follows: 


r" sin(0" - i0) 
sin 0/ 


(34) 


Since subroutine MODATM is called by subroutine REFRAC and since subroutine 
MODATM is called last for the altitude corresponding to the middle of d, the array 
ANS may be used externally to calculate the amount of water vapor or the total atmos- 
pheric mass that was traversed over distance d. For the initial calculation at the 
target point, the angle £(i.e., 0”) is needed; therefore, subroutine PATH is pro- 

vided to calculate £ for the programer. 


Subroutine PATH 

The principal purpose of subroutine PATH is to calculate the angle however, 
while calculating £, it is also convenient to calculate the columnar mass and the pre- 
eipitable water vapor along this path. These three quantities are stored in the array 
ANS. If subroutine PATH is called prior to the calling of subroutine MODATM, ANS(l) 
will be set equal to -1.0, and subroutine MODATM will be called such that subroutine 
INPUT is activated, eliminating the future need to call subroutine INPUT externally. 
Subroutine PATH is thus programed to be called only once for each radiosonde 
sounding. 

The initial guess at £ is calculated by finding (r ' - r,'), the vector from the 

sp 1 

target (1) to the spacecraft (sp), as shown in figure 2 and as developed in reference 10. 

The components of (r ' - r,') are 

sp 1 


Kx = ( R e + Z sp) cos 9 sp cos % - ( R e + Z l) cos 9 1 cos *1 <35> 

R Y = ( R e + Z sp) sln 9 sp cos % - ( R e + Z l) sin 9 1 cos <36) 


and 


R Z = ( R e + Z sp) sta * S p - ( R e - Z l) sin * 1 < 37) 


14 



where e is the longitude of the spacecraft, 0 ^ 

is the latitude of the target, </> is the latitude of the spacecraft, Z 

sp 


is the longitude of the target, <£. 

is the altitude 


of a spacecraft above the earth, and Z, is the altitude of the target above the earth. 


The components (R^, R_, and R^J are found by coordinate transformation in 

the coordinate system of the target to be £", tj", and which are the respective 
distances southward, eastward, and upward from a local station to the target. 


r 

7] n 

C” 


'ib 

sin <£, cos 0 , sin 0, sin -cos 4>, 


1 

-sin e 


1 


1 

COS 0 


1 


cos cos cos sin e ^ 


0 

sin <p 



Kx 


r y 


R z_ 


(38) 


The unrefracted zenith angle 


£ = tan 


-1 


V(r ) 2 + W ') 2 


C” ■ 


can then be found. Next, the angle TOS between r ' and r g ' is calculated by using 
the definition of the dot product " 


TOS = cos 


{ 

r ’ \ 
sp 

\ r 1 
\ sp 

r l / 


(40) 


so that the best refracted path from the target to the spacecraft (fig. 1) may be found 
by iteration. 

Iteration of paths from the equations developed in the description of subroutine 
REF RAC is used to find 0' and i j/ for each level, and since 

A£ = <P" - (41) 

integration proceeds until 


2 AZ = Z - Z. (42) 

Sp A 
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Then. 2 A£ is compared to TQS for the purpose of iterating on £ as follows 


?(t + At) = £(t) - 


(2 Ag - TQS) 
2 


(43) 


until js A| -- TOSj 0. 0001 radian (0. 0057°). This procedure yields an accuracy on 

“3 ' 

£ of approximately 3 x 10 radian (0. 17°). The quantities columnar mass and pre- 

cipitable centimeters of water along this refracted path are calculated, respectively, 

in the following equations. 



d 


(44) 


and 



1 


qp ds ~ 2 qpd 


(45) 


The increments on AZ are made to be multiples of 10 smaller than Z 
that 


sp 


- Zj, such 


Z - Z. = AZ • i 
sp 1 


(46) 


where i is 10, 100, 1000, et cetera and AZ < 0. 2 kilometer. 


Subroutine ATM0S3 

The subroutine ATMGS3 reproduces the U.S. Standard Atmosphere, 1962 
(ref. 1). Subroutine ATMOS 3 is called with geometric altitude from which geopotential 
altitude is calculated. The equations which are subsequently used for ATMOS3 are 
many of those developed for subroutine MODATM. Equations (3) to (5) and (8) to (15) 
are common to both subroutines. The main difference between subroutines ATMOS3 
and MODATM is that in subroutine ATMOS3, all the significant levels are included in 
a data statement so that no data cards are necessary, and the output variables are 
more limited; that is, only the first eight variables in array ANS are available. 
These variables are pressure, temperature, density, speed of sound, acceleration 
of gravity, molecular scale temperature, molecular weight, and coefficient of 
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viscosity. The main purpose for including subroutine ATMOS3 is that if atmospheric 
data above the maximum- altitude radiosonde data are required of subroutine MODATM, 
then ATMOS3 is automatically called. The main impact subroutine ATMOS3 has on 
analyses is that if the maximum usable radiosonde altitude is <10 kilometers, signifi- 
cant water vapor will be ignored since the subroutine ATMGS3 includes no water vapor. 
Instructions on the use of subroutine ATMOS3 are included in comment cards in the 
subprogram. The computer print-out, including all subroutines, is shown in 
appendix B. 


CONCLUDING REMARKS 


It is hoped that this nonpredictive model atmosphere for earth resources applica- 
tions will fill the need for atmospheric data until predictive postflight or flight models 
can be developed. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, November 15, 1969 
160-75-03-00-72 
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TABLE I. - LAKE CHARLES, LOUISIANA, RADIOSONDE AND CODE 


May 10 1969 0000Z 

TT 60004 72240 99016 23266 01008 00146 21467 00512 85517 
08463 35017 70118 04273 32033 50577 13571 29543 40743 2656S 
27572 30946 38567 27590 20217 519// 15400 589// 10650 673// 

88999 

66280 275950 

VV 6000/ 72240 00016 23266 11970 18068 22831 06662 33813 
11075 44609 02171 55400 26569 66290 40166 77243 461// 88227 
451// 99193 535// 11100 673// 31313 25069 451// /////0 

QQ 60000 72240 90012 01008 35512 35007 90346 36009 36013 34524 
90789 33530 34031 33031 91246 31535 32539 31534 9205/ 29044 
27582 9302/ 27588 275950 


2nd Trans 


WW 6000/ 72240 70866 661// 50071 633// 30391 551// 20653 
497// 10115 411// 07358 403// 

88950 681// ///// 

779990 

YY 6000/ 72240 11950 681// 22920 657// 33600 665// 44230 
511// 55100 411// 66070 4O3//0 

LL 60000 72240 XMTD0 


a The significant level code is VV. For VV, the code is 
iippp TTTdd where 

ii = identifier of a set of data; the two characters are 
identical (e. g. , 00,11,22,33). 
ppp = pressure in mbar except the 4th character from the 
right is suppressed (e. g. , 970 = 970 mbar, and 
016 = 1016 mbar). 

TTT = temperature, + if last digit is even, and - if last digit 
is odd. 

dd = dewpoint temperature. If 00-49, multiply by 0. 1 for 
°C; 50 = 5.0° C; 51-55, not used; 56-99, subtract 
50 for °C. 

That is, 02 = 0. 2, 56 = 6. 0, 60 = 10. ) 

Slashes indicate no data. 


TABLE n. - INPUT DATA 




CONTINUATION! 









Figure 1. - Refraction-path geometry through a spherically 
symmetric atmosphere. 


z 



Figure 2. - Resultant vector from the target to the spacecraft in fixed- earth 

center coordinates. 
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APPENDIX A 


DERIVATION OF VIRTUAL TEMPERATURE T : 


The equations of the state of dry air 


( P - f w e ) m d 


p d “ RT 


of water vapor 


f em 
w w 

P w " RT 


and of wet air 


w w 
p “ p d + p w “ RT + 


£ em_ (P - f w e)m d 


RT 


can be used with the mass percentage formula for molecular weight 


m = 


100 


M i p w + p d 


m. m m , 
1 w d 


to give a formula for the relationship of temperature, molecular weight, 
water-vapor pressure 


f em + (P - f e)m, 
w w \ w / d 


m = 


RT 


f e + /P - f e\ 
w \ w / 

RT 


(Al) 


(A2) 


(A3) 


(A4) 


pressure, and 


(A5) 


24 



Equation. (A5), when simplified, becomes 



1 - 0.37803 



(A6) 


By employing the definition of T* 


T* 


m ,T 
d 

m 


(A7) 


and by using equation (A6), the exact expression for T* may be found in terms of T, 
e, and P 
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APPENDIX B 


SUBROUTINES 


QULQJ L»_ SUBROUTINE MOUA III < / . PP , T EST . X L A M U A > 

□ 0103 2* DIMENSION H « 25 ) » P < 25 ) , T ( 25 } , To < 25 ) , AnS ( 35 1 » T V { 25 1 

00104 3* _■ COMMON ANS . _ 


OOiOS 
QQ1QS 
OOiOS 
-M 1 .0 5 


4 • 
5® 
6 ® 
7* 


DATA RQ/8.3i432E+U//,XM0/2S„9664/,BETA/l . 4 S3 E-06 / , S / I I Q , 4 / , RE / 6 , 37 
LL£9?E *03/ ,G/9SQ«665/ , CONN/ s 3 ,41 6 3 1 ? 1 7E-Q2/ 


00 105 
QQ 105 


«® 

9* 


OOIOS 

OOIOS 

10* 
1 1® 

Tfo 1 05 

12* 

00105 

1 3® 

OOIOS 

14* 

00 105 

IS® 

00105 

1 e»® 

OOIOS 

17* 

00105 

1 8* 

OOIOS 

I 9 • 

OOIOS 

20® 

OQIOS 

21* 

OOIOS 

22* 

- OOiOS 

2-i* 

OOIOS 

24* 

OOIOS 

25* 

00105 

26 ® 

Q0105 

27* 

00105 

2B* 

OOIOS 

29* 

00105 

30® 

OOIOS 

31* 

OOIOS 

32* 


C 

C 2 IS I N K M. PP IS IN MB ' 

C ANS IS OUTPUT VARIABLES 

C XLA MQa. IS THE WAVE LENG TH .IN MI CKO NS FOK ftHICH YOU aR£ CALCULATING 


ATM CIPHER IC REFRACT Ion 

IF TEST ,EQ. pres THEN PRESSURE IS USED AS HEIGHT IN 0 1 C ATOR 
IF TEs't • NE • PRES THEN GEOMETRIC ALTITUDE { KH ) I HEIGHT i NDTcATOR” 

YOU MUST S t T ANS 1 I I *• 1*0 0 £K ORe £ NT £ RINg THE SUBROuT|N e the FIRST Time 
RU IS THE UNIVERSAL GaS CONSTANT BASED ON THE CARBON j 2 aTOMjC HEIGHT 
lcAk£_ JJL_tR 5S / 4 OE G .. R EL V IN. 7JS « - H&LEj 


XMO IS MOLECULAR- WEIGHT OF A l R CALCULATED FROM THE COMPOSITION OF D«V 
AIR USING THE CARBON 12 ATOMIC WEIGHT SCALE. FUUN 0 IN ThE U. S. 

STANOARD ATMOSPHERE 1962, PAGE 9» GIVEN IN GM/TgM-MOLE) " 

beta is A CONSTANT used IN SUTHERLAND’S VISCOSI TY EQ UATION. GIVEN IN 
KG/SEC-M-IDEG KElVIN**i/2) 

S IS SuTH ERL aNP*S CONSTANT^ In DEG, KELVIN 

re * the mean Radius uf the earth in meters as given by the smithsonian 

METEOROLOGICAL T A BLtS, S|XTH EDITION, PUBL I CAT I Qn_4UI 4 , g. J, 

LIST, 1966 

G IS ACCELERATION OF GRAVITY A? 0 £«U I POTENTIAL SURFACE LEVEjL G I VEN IN 

CM/SEC**2 

CONN js A CONSTANT GIVEN A S -M*G/RQ WHERE M IS MASS aNd 6 aNd RO ARE AS ABOVE 



aoios 33« c the foulowing is an example of a calling program for mqd a th an d path 

OOIOS 39® C DIMENSION ANS135) 


OOiOS 

36 ® 

c 

XL AMD A* • 6 

00105 

3 7 * 

c 

2S«20«U 

OOiOS 

3U® 

c 

PHlS»30oQ 

00105 

39® 

c 

THETAS-9Q.U 

00105 

< 10 ® 

C 

ZL»0«0 

OOIOS 

_JLLl_ 

c 

_ P.Hlk-^0.0 _ 


OQ 1 OS 
00105 


92» 

* 53 ® 


T HE V AL“?0 » 0 

CALL PATH (XLAMtM, *S, PHIS, THETAS, ZL,PHlL,TH£TAEl 


00105 
QOIQS 
00105 
00 1 ns 

Hi® 

95® 

96* 

97* 

c 

c 

c 

c_ 

WRITE (6,3) ( ANS ( M >K>2 1 ,23 ) 

3 FORMAT ilXi///»lX,lP3E!H, < 4l 
TEST-1HPRES 
Dn 1 1=1 ,20 

OOIOS 

98* 

c 

Z*l 

QOIQS 

i»9* 

c 

PP«1000<!Q-Z*50, 

QOIQS 

50* 

c 

call modatm (Z.pp.test.xlamdai 


QQ LOS Si* 
OOiOS S2» 
QOIQS S3* 


Q010S 

S9 # 

OOIOS 

55* 

00105 

56® 

00105 

5 7 * 

00115 

58® 

QO 1 1 S 

59 « 

00156 

6 0® 

00 1 20 

6 i* 

0012 1 

62® 

00122 

6 3® 

QQ129 

69® 

QQ129 

6 5* 

0012S ■ 

6 6® 

00130 

67* 

00131 

68 * 

0013 A 

6?* 

Q Q 1 3 6 

7U * 

00137 

7 1 * 

001 AQ 

7 2® 

00191 

73® 

QQ 1 9 2 

79« 

00193 

7 5® 


c 1 write < 6,2 > test » Z • Vans < n f , n* l , 29 ) 
c 2 FORMAT ( l X » A 9 , A X , 1 P i 2E9 « 3 , / , J P 1 3 E9 . 3 » / / ) 
_£ C ALL e xit _ ' __1 _ . 

C END 

c 


c . 

CT*28e» 15/273. 16+1 * U 

C cr IS \ . Q * RATIO OF SURFACE TEMPER ATURE TO ICE T£ HP£RaTURe 
I F ( AN S ( I ) . GE . U * 0 i GO TO 1 5 

ANS ( n*o.o 

call input (p, t»td,h,tv,m) 

IS IF <TEST,EQ«9HPR£S) go to 7 ... 

HA® RE*2/tRE+2)®lOUU»0 

C HA IS GEOPOTENT I AL ALTlT U DE IN METERS 

23 DO II 1*1 ,M 
I I * I 

1FIHU! “HA ) 11,12.13 

1 1 CONT INUE 

9 CALL ATMQS3U) 

A N S ( 9 ) = 0 » 0 ; 

GO TO 52 

13 1 = 11-1 

DH = H 

0= 1 TV (l* l 1 "T V 1 I > * /OH 



001*4 4 
00 145 
00146 

00146 

00146 

00 14 6 
0014 6 


7b* 
7 /• 


78® 

. 2 . 9 * 

au* 


a 1 » 

82* 


W>| f { 1 + 1 ) ” T ( I ) I / Utl 
DW»(T0(1+1)-T0(1i)/0h 
0 H xH ( 1 1 -HA 


c 

c ® * < 

C 

c HEIGHT 


» H • ' IS"Ttr«£T£«S 


00146 83® C HEIGHT »Z* IS IN KM 

00146 84* C 

00146 8S*_ C ANSI l ) IS. PRESSURE _ 

00146 86* t PRESSURE IS IN MB 

00146 6/* C 

00146 88® C ANSI 2) IS TEMPERATURE 

00146 81® C T EMPERATURE JLS IN DtG KELVIN - 

00146 90® C 

00146 11*. C ANSI 31 IS DENSITY _ J - — - 

00146 92® C DENSITY IS IN SM/CC 

00 146 g j« C — 

00146 94® C ANSI 4) IS SPEED Of SOUND 

00 1 46 . 95* JC SPEED OF SOUND IS IN 11/SeC 

00146 96® C 

002 46 97® t_ANS< 51 IS ACCEEERaIION QF GRAVITY - - - 

00146 98® C ACCELERATION OF GRAVITY IS IN C«/S£C®®2 

00146 99* C 

00146 100* C ANSI 6) IS VIRTUAL TEMPERATURE 

001 46 SOI® C TEMPERATURE IS IN DEG kElVIN ... . .. 

00146 102® C 

0Q146 JOi* ._£ ANSI 71 IS MaLE.CULAK HEIGHT — 

8B1U 181® i ansi 8) is coefficient of visc osity . 

00 146 TOfe® rVTsCOSlTY IS IN KG /(M SgCl 
00146 i 08* C ANSI 91 IS DE» POINT TEMPERATURE 

Q014& 1 0 It C TEMPERATURE IS. IN DEG KELVIN 

00146 110* C 

00146 til® C AN SI 1U1 IS . M] X IN G RA TI O R 

00146 1 12* C MIXING RATIO IS IN P aR T S / T HOUS AND WE. (0/001 GM/K(i 

0.0. 1 46 _I13®_ C... , - - - 

00146 114® C ANSlill IS SATURATION MIXING RATIO RS 

0 014 6 J IS* C. SATURATION MIKING RATIO IS . I M._£ AR I S/.T HQ U S ANlQl I*JEjl_ 10/0.0* .SM/JCG. 

00146 116® C 

00146 1 1 7 ® C ANSI 12} IS RELATIVE HUHIplTY 



001*4 6 11*3* C RELATIVE HUMIDITY IS IN PERCENT *,0/0j 

00 I *46 11?* C_ . __ _ .... ■ : — 

GO I <16 120* C ANSI 13! I 5 SPECIFIC HUMIDITY 

0(11*46 121* _ C SPECIFIC HUMIDITY |S JN Ofl/KG ' 

00 1 R6" “ 122® C ' 

001*16 123* C ANSUHJ IS SATURATION SPECIFIC HUMIDITY 

00 I <16 1 2*4® C SATURATION 5FFCTFTC H‘JM i o TTY Ts In 5m /k G 

OOIH C — 

001*46 126® C ANSI IS) IS PRESSURE SCALE HEIGHT 

00 I <46 12 7® c PRESSURE SCA L E HEI GHT IS in KM 

00 I *4 6 128® C .• 

0QH6 129® C ANSI Hi IS DENSITY SCALE HEIGHT .... 

001*46 130* C DENSITY SCALE HEIGHT IS JN KM 

0 Q1 H A 132® C ANSI 17) IS REFRACT 1 VE INDEX DEVELOPED UY EDLEN IN ftWs OF y^WEXENGTH ALONE 

001*46 133* f INDE X i S F OR aIR a T 2 8 8 K£LV jn AND 7<j U.MM HG __ 

Qqj*4& 13*4® C 

001*46 1 3 S ® C ANSI 18) IS REFRACTIVE INDEX DEVELOPED 0Y PLNNDORF JN TE«MS Of 

001*46 136® C WAVELENGTH » TEMPERATURE* AND PRESSURE 

001*46 1 37® C ' _ : 

001*46 136® C ANSI 19) IS THE WATER V A P 0 R PRESSURE I N lit* 

001*46 139® C .• i 

001*4 6 r*4U® ' C ANST20i I S ThE S aTuRaT i UN WATER V APOR PRESSURE IN MB 

00 1 *46 i *11® C _ _ , 

001*86 1*42* C ANS 1 2 1 { IS THE 2£N|Th ANqeE FROM SROUNOSTAT ION IN. RADI aNS 

00 1*46 1*43® C _ _ _ . 

001*46 1*4*4® C ANS 122 i * THE TCiTAL GM/Cm»»20r COLUMNAR MASS ALONG THE SlaNT^FaTH® 

00 1*4 6 I «4S® __C 

00 1 *46 1 *46® C ANS123) * TOT AL . GM/CH * * 2 OF WATER VAPOR ALONG THE SLANT PATH* I? I.S 

001*46 ..... s *4 7 *_ c equivalent to pRtc.ufi table cm of water 

QQ1*S& 1*40® C 

001*46 J*)?® C ANSUSJ » IQTAL PATH LfcNfiTH IN CM .... — . . 

00 1 *46 ISO® c 

004 446 IS L* . C A 14 S I 2.1.1 THRU ANS I 2* 4 ) A RE C A L C U L A T E P . IN... SUBR QUTINE PATH . 

00 8 *46 15 2® C 


001 *46 153® 

CO I *46 15*4® 

QO i *4 7 1 5b * 

00150 156« 

00.151 1 S7» 

150* 



i 0 $> $ 0 0 $ $ 0 0 0 ft ® ® 0 $ $ 0 ® 0 ® ft ^ 0 9 0 0 # 0 # w # f.' # 0 0 © 0 0 0 


ANSS2)«TU.l"l8»0H 
ANSI 6 4 *TV < I 1 _ U«Dh 

ANS4 9)*»TPi 1 ) -PW*DH : 

ANS m»PRES(P u I ,D ,i V ( l > ,ANS (6 | ,-DH ) 


159* GO TO 1*4 


00152 
001 53 



C.' C i & H 


1 60* 


00 155 

00156 

00157 
00 1 60 
00161 
00162 
00 163 
00164 
00 1 65 
00166 
00467 

00170 

00171 

00172 

JULLU. 


QQA74 

00175 

00176 

00177 
00177 

002 Q 1 
00201 
00202 

00203 

00204 

00204 

00205 

00206 
00207 
00207 
00212 


161 * 
16 2 ® 
l 6 3 »< 


1 6 4 * 
1 65* 
166* 
1 6 7 * 
166* 
i 6 9 * 
170* 

m* 

172* 

173* 

174* 

-i 75* 


176® 
17 7 • 
178* 
17?® 
180* 
181 * 
182* 
183® 
184® 
188* 
186* 
187® 
188® 
18V* 
1 9 U • 
191* 


00213 

MOZ 16 
00220 

00221_ 

00222 

002 25 

00226 
00227 
00231 
00233 
00235 


1 92® 

. 1-9-3® 
194* 
,195* 
1 96* 
19 7* 


12 


-1 4 


82 


9 0 


1 * I 1 

AilS < 1 1 »P ( 1 ) 

A N S ( 2 ) * T ( 1 1 
A N S ( 6 1 sT V (I ) 

A N S 1 9 1 * T 0 l 1 , > 

ANSI 5 > *G* ( RL/ ( RE + Z 1 > **2 
ANS131=ANS( 1.) • xMO / ( Ho • ANS C 6 1 ) * 1 000 ■ 

ANS ( 4 ) sSQRT 1 1 . 4 *RO * ANS 1 6 ) ■/ XMO ) / 100 , 

ANS ( 7 1 =RMO»ANS t 2 ) / ANS t 6') 

ANSI 8 1 *8 E I A * ( SQR T_( ANS ( 21 ) 1 **3/ 4 A NS ( 2 ) * S 1 
ANSI 1 9 1 "E l ANS l V > f 

ANS(20)«E(ANS(2> ) 

ANSI 1 0 >»R < ANS (l? ) . ANS < 1 1 ,ANS(2 ) > 

ANS U l 1 *R < ANS l 20 ) , ANS U 1 t ANS < 2 1 1 
ANSI 12)*{ANS( 1 0 1 / ANS ( 11 1 1*100.0 

AN SI 13 1»Q( A NSI i > , AN5{ 9) ) . 

>Q ( ANS ( 1 ) , ANS ( 21 ) 


ANSI 141 , , 

ANS U5)-RQ*ANS(6»/(A Me *ANS(S> 1 * l.OfOb 
ANS* 16) *ANSU5) / 11 »0*R0/ 1 XH0*ANS(5| 1*0* *0 1 1 

IF (xLAMOA.GE. 1.2500^001 GO TO 30 

THIS MEANS IF XLaMqa IS .GE. 1.25 CM USE MICROWAVE KEFRACTlVjTY 


, n t, M II *> 4 r ^ t- AA » ' w i mf e Uk. * w j *- II S V'V" n » i. -’»«v 

AN9 ( 17)» L» y* 1 • QE-0.8® (643-2.8*294 ?J_1 0 * / 1 1 4 A * - 1 * / (XlA m QA*» 21 l_*£SS4CL*Z 

1 (41 .-l ./( ALAM0A**2> 1 1 

ANSI 18 1*1 .0*1 ANSI 1 7 1-J.01 MOT/ ( 1.0* AnS (21/273. 16} 1 *ANS ( 11/1013.25 _ 
GO TO 31 

30 ANSI 111 *1 .C+l .0E-06* (7 7.6*ANS( 1 1/ 1 ANS (2 1 1 ♦ 37 30(10 . n* ANS ( f 9 1 / l.ANS ( 2 

1 1 * » 2 ) J “ T " 


198® 
1 99 * 
200 ® 
201* 
202 * 


ANSI 1 7 1 * A N S ( 18) 

31 RETURN 
7 00 16 I - i , M 
PRESSURE 

1 1*1 

IF (PP-P( I 1 ) 16 ,4 i , I 7 

16 CONTINUE 
HAaO.O 

OHA* 100*0 

5 1 00 4 8 I » Is 11 

H A«H A*OH A 

CALL ATM0531HA) 

IF (ANS( 1 l.LE.O.O) G 0 TO 42 
IF ( ANS ( 1 1 . LT . PP1 GO TO 49 

48 continue 
42 00 10 I ®2 s 35 


Q02JiQ 

2 n 3 * 

IF 1 I .EQ.21 .OR. I.EW22.0R.J 

i t E.Q « 2 3 . Q R «_.I * E S «. 291 GO TO la 

00242 

204® 

ANSI 1 1=0.0 



00243 


20 b * 


» 0 COWT 1 NUE 



Q02H5 206* ANSI I 7 > * J * 0 

CJ Q 2 *4 6 .. . 2U.Z* — _ AMS llfl I * i • Q — - 

002*47 20&® 2»«A 

00250 209# RETURN _ _• 1 

00251 210® ‘M I I >*RE/ ( 1000.0* iRE’HU I > /lOOO.O ) > 

.00252 2l!i. . ' SaJO 12 - 

00253 2J2® *49 IF S ABS ( ANS ( 1 > -PP > *Lfc* ( .00i*PP) • So TO 50 

QQ255 . .. 2 i i*L HA«.H A<*DHA._ - 

00256 2 1 *4 * DHA“DHA/1Q.Q 

00257 24Jsji _fiO — LQ-JLJ — - — 

00260 216* 50 2®hA 

. 0.0.261 24.1* GO TO 9 

00262 2{U« 17 I »i 1-1 

.00263 2l»* ; tLMXlil tlL-LSLl U .. - 

0026*4 220* IF(O) 20,21 i20 

00267 221* 20 0 «c 0 NN/At 06 (P ( I»M /P( in»AtQG(TV(l»n/TV(lt) 

00270 222* ANS<6»«TV< U*(PP/P< 1 1)»*<0/C0NNJ 

0027 1 223* HA«H (1 ) *( A NS <6)-TV( 1)1/0 

00272 22*4* GO TO 22 

00 272 225* C HA IS [N METERS 

00273 226® 2i Ha*H ( J ) *T V (I) • ALOG < P p/P4 I )1 /CONN 

00 27*4 22 7* 22 2® HA»RE/ ( 1 OOP . 0* 1 R E" H A/ 1 OOP « Q ) i 

00275 226* GO' TO 23 

00276 . ?y«_. . END ... . . 

END _0f UN I V AC 1 1 08 FORTRA N ’ V COMPILATION, 0 *01 AG NOST IC* MESSAGE ( g ) 





t<S 

!SS 


00105 

1 9 

OOiQS 

2 6 

00101 

3 9 _ 

00 J03 

09 

00 1 00 

& * 

00 s OS 

6 # 

00105 

7® 

00 105 

e« 

00 1 OS 

9 ® 

00105 

10® 

00105 

i l® 

00 1 05 

12* 

00105 

i 3® 

00 5 05 

1 0» 

, 00105 

15® 

00105 

1 fe® 

00105 

1 7» 

00105 

IS# 

00105 

1 9* - 

00105 

20® 

_QP5 05 

■ 21® 

00 105 

22® 

00)13 

23® 

00113 

20# 

00113 

25® 

00113 

26® 

00113 

27* 

008 13 

28* 

00113 

29* 

00113 

30* 

00U3 

31 • 

00113 

32® 

00113 

33® 

00 l 13 

30® 

00113 

35* 

00113 

36* 

00 113 

. 37* 

00113 

30® 

00113 

39# 


SUBROUTINE ATMQS3 (2) . .. 

C SUBROUTINE FOR THE 19 6 2 STANDARD 

C l I5» altitude.. JN KM . 

DIME NS JON hI23)*T(23),P1231»aNS.<3S)»a<23>,ZZ(23» 
COMMON ans 


"DATA H/-50GQ. » □ •a,'llUUQ*Q»2OQOQ*O»32(j£|a*O»97QOll*U»5 2 0UQ*O»6lO®.O* u » 
l 7 ¥ 0 00 o , 58700.2, 90052, ,108129.8,1 1 7 7 7 7 . 7 . 19 65 9 3 * 8 . 1 S&O 7 3 * 6 , J 65S?A»1. 

2, 180088.55,22 1 9/2. 606 ,286086.09 ,37633 1 *361 » 063556*05 ,508275.86 » 

2630590.90/ ,1/320. 65* ' _ T *_ 

3 » 208 o 15 »21"6*6S»216.6S >228«65»27 0.6S»27 0»6 5»25'2»fS , l80»65,ja0.6S, 
9210.65,260.65,360.65,96 0*65? U J 0 . 65, 121Q. 6 5 , 1350.6 5*1 550,6591830.6 
55,2l6Q*65»292b«6S»2&90»65«27Qo»65/,P/1.77687E*03jl.01325£*Q3» 

62»263.20E*02 ,5.07487t®a \ » 8 » 68Q 10,1 « 10905 t 5»90Q Q5£- Pl t i »3 2Q99E-0J , 

7 i*037 7E-02t t . 6 0 38~E-U3 , 3 # QQ75E-Q9 9 7 « 35 9 9£-05 , 2 . 52 4 7g-05 , 5 . 06 I 17E-06 , 
83*699 3E“0 a ,2»7926E”06 . L . 68 52E-Q6 , 6 • 9 609 E-0 7 , l * aH 3 3 E* 0 7 , 9 . U3O0E -Q6 , 


9 l . o 9S7E“08»3.9 502E-09, l« 19 {§£-09/ , A/320.8&0* 

1 286* 15,216* 65,216,65,228 9 Jb 5 t-- 27 0 • 6 5 >27 Q * 85 » 252 » 9 6 , 


' 65 i 1 8Q ,65 , 


2 210.02 ,257.0,399,99 f 8<?2,79 » 1022,2 , 1 105.5,1205,5 , 1 3<! I .7 , i 932 , J , 

3 5 987.9,1999*2, |506.i , 1.507*6/ ;ZZ/ -fcOOC . ,0 • C » 11 QQp* » 20QQ3* ,32000. , 

907000., 52000. ,6 J 000 •',79000,, 9 00 00. »i 00000. »1 10000* ,120000= » 5 50000, 
5, HGQOO* , 170000*, 190000. ,230000* ,300.050. , 0 0 PS Ofi* lSjoMQ CLt ^ - 

6700000./ ‘ 

DATA S./ 1 1 . 0 . • 9 / , C 0 N N / “ 3 * H 1 6 3 1? J 7. E - 0 i / » 8 E / 8 *. 3 6E ±M/ — 


C 22 JS THE GEOMETRIC ALTITUDE FOR BREAKP OINTS A B OVE ^ . 

__rf— — T ~~ A uf IN S£0P8 teNT l A t METERS FOR SIGNIFICANT L fc V £ L 5 

c 0 IS THE TEMPERATURE GRAfl JENT JN THE VERTI CAL (PES/Ofe Op WI 

C TUI IS THE MOLECULAR SC A lE TEMPERaTuRE AT A SIGNIFICANT LEVEL 

C Mj) IS THE KINETIC T EHPeR A.TUR£_j|T. UHlUl 6 NJ FJ£AilI_tfcWEkS — 

c Pdi IS THE PRESSURE IN LfcS/FT**2, ACTUALLY IT WONT MATTER AND PRESSURE CAN 
C BE IN ANY SET OF UNITS SIN CE ONLY THE RATIO AT VARIOUS &L T !T U PES> RELATIVE 

_ f 0 ~"p ji j j s U sed ~ 

C ANSI!) IS THE RATIO OF PRESSURES (P/PSi.1 _ 

C ANSI 11*1 .01 325E*Q3 FOR PRES SM MB 
r A N S 1 2 1 IS THE .RATIO OF LeMP£RAJ.U.R£ (T/TSL) 


c ANS 121*288.15 FOR TEMP I N BEG K 
C AN S 1 3 1 IS THE RAT IO OF O enSITIES 


GO i 1 3 009 C ANS<3»*I.225E-03 f OR DENSITY IN GH/CC 

OQ 5 5 3 tlj C AN St9» IS THE RATIO Of. SPEED Of SO UND tC/CSLL 

OO 543 02» C «N519{#39U,290 FOR SPEED OF' SOUND IN M/Sec 

005 J 3 ......03*! C. ANSI 5 ) IS THE ACCELERATION OF .GRAVITY. . t<»/ 6 .SL J_ 



001 13 

~ Q1LLLL_ 

00113 
. 00113 .. 
00113 
001J3 
004 1 3 
001 13 
001 1 3 
00113 
00113 
00113 
001 1 3 
00117 

00120 

00 . 123 . 

00124 

- 00125. 

00126 


99* 

95* 

96* 

*4 7 * 

9 8 # 

9 9* 

SO* 

SI* 

1 * 2 *" 

53* 

59* 

55* 

56® 

5 7* 

56* 

59* 

60 * 
61 # 

6 2 * 


t A NS I 5 ) * 9 80 s 66 5 FOR ACC °F GRAVITY IN CM/ ( SEC*®2 1 
..RAI11LJ1S... MOk.tiJ»I.LM.-.SLCA.Lfc - 


iMU/MUSLl 


AN5 1 6 ) *288 * 1 5 FOR TEMP *N DEG K 
.ANSI?) IS THE MOLECULAR' »fcl<SHT 
ANS ( 8 ) is THE RATIO OF CQ.Ef Of VIScQSlO 

ANSl8)*».7899E-05 TO COEp IN KM/M-SEC 
» IS THE VERTICAL- KlNETlc TEMPERATURE GRADIENT 
T HI S_Ra 01 u S «RE « IS C H OS£ N Tp A6R.EE 1 IT H THE u 

ALSO fs A BEST F I T 0 ALL "LE VELS 8ELg» 90 KM® 


s standard at ho km, but it 


Asoyt 9o km the .levels 
thit.aRe break Points here calculated fro m geo metr ic to se o p_ u sjlN g _ »_R£». 




Z«z* 1 0 0 0 . 0 

IF (Z-7QQ0Q0.Q 1 1 0 s 5 q 8 50 

10 continue 

Ha«RE*Z/(RE+21 

ANS151»RE«*2/ s (R£*Z »**2 1 

DO I M* I 8 23 ' . 


^RTSTn^l 

63* 





00132 

69* 


IF- SHU 1 "HA ) 1,2,3 



JJQ135 

6b* 

i 

CONTINUE 



00137 

66* 


GO TO 50 



_ .00.190 . . 

- ii z*.... . 

.__A 

...1*1" 1- - - - - - ...... _ - . . . 



00191 

■ 68* 


D»|TI|*ll*TIII)/lHllHl'H(II| ■ 



00 1^2 

69* 





00193 

70* 


60 TO 9 



Q0J99 

7 1 • 


A N S 1 6 J.« T | l 1 / T ( 2 1 . .. _ _ 



001*15 

72* 


ANS ( 2 ) * A l I 1/A12) 



00196 

73* 


D«|T!I®ll-TSil)/!H(Ul!-H{ii) 



00197 

79* 


GO TO 5 



00150 

75* 

9 

IF <90000. 0-Z> 7,7,9 



00153 

/ 6 * 

; 

ANS‘6l«<TUI-(TU*l|-TU)»/lZZ(l + l)”ZZUl)*U2Ul-Z>l/T(2> 



.Q01A9 . 

.77* 


_.ANSJJLUL1 A.1 I J -...LA i L* i. » - A (. I 11/ LZZtl + 1 J.-^iIllJ_*.LZi.U-.L-Z.iiL/AlZl . - 



. 00155 

7 a * 


GO TO 5 



0.01 56 

7 9® 

9 

ANSl6)*(TU>-D*(H(I<-HAn/T(2) 



00157 

BU* 


ANS ( 2 1 * ( a ( l 1 - »« • ( H l I * - H A » ) / A ( 2 1 



0016U 

a i * 

b 

IF ( 900Q0.0-Z 18,6*6 



GQ163 

0 <2 * 

6 

ANS( 7 1*28.9699 



00169 

S3* 


GO TO li 



00165 

3 ^ * 

8 

ANS(7)»28«9699»AnS(2)/ANS(6) 



00 16 6 

85* 

-J-J 

ANS ( 9 1 ® SGRT(AnS(61 * 





03 


0 0 16 7 

6 6® 

QQ 1 70 

87* 

001 73 

86* 

0017*} 

89® 

00 175 

9 0® 

00176 

9 1 * 

0 Q 1 7 7 

92® 

00200 

93* 

00201 

98® 

00202 

95* 

00203 

96® 

0020H 

97® 

00205 

98® 

00206 

99® 

00207 

100® 

00210 

101* 

0021 1 

102® 

00212 

103® 

00215 

1 o*l« 

00217 

IQS* 

00220 

1Q6® 


A N SlB) = ((T(i)<-b>/lANs(2)*T(2!^S»i*SiJHTCUNSi2!) s *J> 

IF I D 1 1 2 , 13 , 1 2 

12 C0NNsD*AL0G(Pl 1*1 J / P ( I > ) / i A L U G I T ( I + i )/T ( I ) i I 
Ai-iS* 1 ) a P I I )/pI2l*<ANs(61®T ( 2 } / X S I I )*®tC0NN/Q) 

go to n 

1_3 COnn=ALOG«PI i-HU )l«j(_II 

ANS I 1 > »P ( i 1 / P l 2 i « tXp l CONN* (ThA-H I 1 M / ( ANS ( 6 > * fTT) ) 1 ) 

lHANSt3)*ANSU)/ANSl6> ___ 

ANSI ()*ANStJ)*l«01325E+03 

ANSI2 >»ANSU)*288. lb 

ANS*3)*ANSl3i*l«22bk.-u3 

An £ [H 1 * An S L 4 !!? 3 ho • 2 ’ H _ 

ANs‘<sT*AN"S(5}®98Q,6fS 

ANS 16 J »ANS (6 ) *288. lb _ 

ANSI8)*ANS(8»*1 .78VA E -05 

Z«2/ 1UU0.U 

GO TO S3 

bO DO & 1 1*1 >8 

5 l ANSI I >®0.0 

b 3 RETURN _ 

END 


end of univac nob fortran v compilation. 


o .Diagnostic* mess a ge{si 









os 

os 


00 1 45 
00147 


00151 

00152 

00153 

00154 
00156 
00156 
00156 
00156 

00156 
00 156 

00157 
00160 


H 0 * 

4 7 * 

4 6* 
iJV* 

50® 

5 l* 
5 2* 
53 * 

54*l 

55* 

56* 

57* 

58* 


1 


IF ( T U ( I 1 * G E • 56 

IF (TQl l 1 I • LE . « 0 1 ) 

TO ( 1 1 *T SI 1 -TD< 1 1 
T ( I ) = T 1 I 1*273.1 6 
TU( I )*TD( 1 1*273* 16 
COnT i n u e 
GO TO 2 


0 .AND. Toil! *LE. 
JO ( l 1 » Tj 1 1 * 273. 1 6 


99,0} T0lIJ*T0U»«-50, 



_C ; 

c 

C T H 1 5 SECTION INPUTS N0N~(.0DED DATA 

' C . _ " _ 

li M»0 

DO 12 1 ° l i 25 


00160 59* C THIS IS THE FORMAT' FOR REAoYN6 "SIGNIFICANT" LEVELS IN nON”CO0EP FORM 

00163 60* READ 15,13) H (I > ,P ( * ) » TU > ,TP 1H _ 

00171 61* 13 FORMAT 1E9.3 ,E1 2. 6 ,(■ 7.2 9.F 3e0 } 

0Q17I 62® C TUill HERE, IS R EL Aj I V E'hOM I D I T V UNTIL. A TO ( I 1 IS FOUNp B T ITER A T ION 

00172 63® . DEL T» 1 00 * 0 

00173 64* GUESS-50.0 


00174 

6a® 

00175 

66* 

00200 

& 7 » 

00201 

68® 

00202 

69© 

00203 

70* 

00206 

71 * 

00210 

...72*. 

0021 1 

73* 

00.212 

74® 

00213 ' 

75® 

00215 

76* 

00216 

77* 

00220 

78© 

00221 

79© 

00221 

ao* 

0022 1 

81* 

0022) 

82® 

00223 

83* 

00226 

8 4* 

00230 

85* 

0023 1 

86* 

00232 

87* 

00233 

8B« 


992 


R 1 »R ( fc"( T I 1 ) » .PI l I * T < l > » 

DO 990 L»1 . 11 ■ _ 

guess=guess*delt 

REL“R1E«OOESS1 (Pill » GUESS I ♦ 1 QO « 0 /R 1 
0«REL*TD( I ) 

IF IQ) 990.99 1 ,995 ___ ' 


990 CONTINUE 
CALL EXIT 

995 GUESS«GUESS-DELT 
DELTsDELT/10.0 

9 9 j IF I ABS ( Q 1 • 6T • • 0 1 I Go TO 992 
T0( I )*GUESS 
IF (P ( I 1 .LE.0.0 1 GO 
M*M+1 

12 CONTINUE 


TO 2 


£*»**•• 

c 

2 


DO 5 l=J ,M 

IFUDtn »LE. Q » 0 1 GO TO 7 . 

TV U 1*T« I 1 / ( i.O“lO.37aO3*E(T0U > )*F <PU > #TU M/PU ) > » 
GO TO 5 * 


TV ( I >»T 1 I 1 
CONTINUE 



GQ235 
. OOiHO 
0024 i 
0024 3 

8V« 
9 0 e 
¥ 1 • 
92® 

2 6 

00 26 I*Ms 2b 

H 1 I J»Ht« )_ - - . 

P ( I 1 ®p 1 M ) 

00 6 I s 1 (H 

00246 
. 00350 

00251 

J9J2252 

00253 

00265 

9 3 $ 
9 ^ ® 


IF i ABS ( T4 11-TD5 i > 1 «6T • 1 #0 + H ( 1 ) * .000/77 > 60 TO 27 
C0NDE»5HCQNI)E ' ‘ - - 

9b® 
9 6 * 


nsati*shnsati 

ON® 5 HON - - - 

9 7® 
9 d » 

27 

24 

WRITE (6,241 H(»)iP( 1»»T(1>, 10(11 «TV(1J , CONDEMNS AT i ,0N 
FORMAT (26X , lP2fcl3.3.QP3F13.2, lX f 3A5l 

002 6 6 
0Q27P 

9 9 * 

J ou« 


IF (CONDE«£0«5H ) 60 TO 6 

C0nDE®5H ... — 

002 71 

00272 

00273 
00275 

101* 

1 Q2® _ 

103® 

104* 

6 

NSATI*5H 

PN»5M - - - 

continue 
WRITE (6,661 

00277 

00300 

10b* 
l 0 6 a 

36 

FORMAT (//l 
RETURN 

0030 I 

107* 


ENO 


ENO OF 

UN I V At 

iios fortran v compilation. o *ujagnostic* message (S) 



03 

OS 


QOIOI 
. 00 1 03 
0010*4 
00 4 05 
00105 
,00105 
00105 
00105 
00105 


00105 

.00105 


00105 

QQ1Q5 


QQ105 
00 1 OS 
00105 
0 0105 
00105 

00107 
00110 
001 1 1 
00112 


00113 
001 19 

00115 

00116 

00117 
00120 
00121 
001.22 
JJU 2.3 
0012*4 


1 * 
2 • 
3 m 
M © 

S « 
6 * 
7 «■ 


subroutine refkal c/i ,Z2,xlamda»phi »phipr»psi , slant j 

dimension A N S 4 3 5 ) 

common ans _ ___ 

data re/63 7 1.299/ 





6 0 
_ 9 _*_ 

1 0 • 

1 1 * 

12 * 

13 * 

1*4* 

15 * 

1 6 * 

1 7 *. 
18 * 

4 9 # 
20 * 


c- in'" order to calculate a continuous path you must extern ally set phT*phjpr 

C 21 , Z j , PHI > AND XLAMDA AKE INP U T VAR IABLES 

C 21 AND 22 ARE IN KM AND xLAMDA IS IN MICRONS , 

C PMIPR, PSI, AND SLANT A«E OUTPUT VARIABLES 

c phj phipri and psi are in Radians and slant js in cm 

C IF YOU WANT AMOUNT OF GM/CM**2 <COL u MNaR NASSl OF ATMOSPHERE FROM 21 TO Z2 
C USE ANS13)*SLANT. GM/CM**2 OF iftATtR IS ANS 4 3 ) *SL aNT »aNS ( 1 3 ) / 1 OOQ ,0, 

C s IN CE AL ’L A NS a RKAV J S I N COMMON , Y&U CA N DO THIS EXTERNALLY. 



C 


S 1 sRE+Z 1 
S2*K£+Z2 


21* __ D£kl“<22-Zl 1/2.0 __ 

22* CALL H0DATM(Z2+DELf .PP.HHALT I , XLAM dA 1 

23* D2*ANS(3) 

2*4* XN2»ANS(1B1 

2b* CALL M p 0 A T M 4 2 1 + U L LT , P.Px,*lliA,LII_»,ALAtl P AJL . 

26* D 1 "ANS 4 3 ) 

2Z*. .. AN 1" A NS 8 18 1 

28* PS1«51NINV<S1*SIN4PHI )/S2) 

29* PHlPR«SlNINV4Sl»SINtPhI)*XNl/4S2»XN 2> ) 

30* SLANT«Sl*SIN4PHl-PSl 1 /SIN (PS I ) * 1 « UE*QS 

3i« return ' ' 

32* ENd 


END OP UN 1 V AC 1108 FORTRAN V COMPILATION. 0 -DIAGNOSTIC® MESsAGElSl 



OOlOi 

1 • 

00103 

2® 

OOIOS 

3 ® 

00 l OS 

•j® 

OQ 105 

5 * 

00 105 

6 ® 

00105 

7* 

00(05 

6 • 

00105 

9 • 

00105 

10® 

00 J05 

11® 

0Q1Q5 

12® 

00105 

13® 

00 105 

1 6® 

OOIOS 

1 5® 

001 05 

16® 


00105 

00105 

00105 


17 ® 
i a® 

19 ® 


C 

C* 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBKUUT I Nfe. PATH IXLAmOA 9 ZS,PHIS»THE! AS*ZL»PH1l,1HEtAL>~ 
DIMENSION ANSiib) ,Alj,3) s B C 3 i »C(3i 
COMMON ANS 

DATA PI/3«i‘4l592dS/»cON/,>0l7 t lS329 2S/,,RE/637Io299/ 


OOaNT 1 T 1 ES ENDING IN S ARE FOR ThE SATELLITE 

quantities ending in l are for the ground local 

-m- AND -02" ARE DUMMY VARIABLES _ 

-AS c YS» AND HS- AKt THE RECTANGULAR coordinates UF the spacecraft 

-AL . V L s AND HL- are THE RECTANGULAR COORDINATES OF l H E GK.OU.NO COCAE 

THE ANGLE ABD IS THE ANGlE BETWEEN THE SUBSATELLITE POINT aND TARGET. 
ANGLE ABD IS found by DSjng The DOT PRODUCT and TAKING JhE INVERSE c os 
• 0092833 RADIANS IS ThE'toTAL REFRACTION ON a P ASS "ThRU U , S . STANDARD 

» su m » | S t he t otal angle change dur ing r e fra ction 

'SUM 1 » is THE SUM OF aLL DELT A ' X I C ALCUL a TED 8 Y La* OF SlNgS 


'SUM2» IS PRECIPITABLL cm of WATER o R G M / CM • ® 2 . 0 F WATER VA POR 
» S U M 3 » is THE total COLUMNAR MASS IN THE SLANT PATH 


00105 

20® 

C 

* SUMR • l: 

ootos 

2 i • 

C 

PHI IS J 

00105 

22® 

C 


00105 

23® 

c 

ANSI2I ) 

0P1OS 

2 **» 

c 


00105 

25® 

c 

ANS I 22 J 

OQ i 05 

26* 

c 


00105 

2 7 ® 

c 

ANS 123 ) 

00105 

28® 

c 

EQU I 

00105 

29® 

c 


OOJOS 

30® 

c 

ANS l 2a j 

00105 

3 J * 

c 


00105 

32® 


00105 

33® 

c 


001 } 1 

■ 3*4® 


PHIS 


aNS i 2 1 ) IS The ZENITH angle FROM gROUNDSTAT ION in KADIaNS 

THE TOTAL GM/Cm«« 2 OR COLUMNAR MASS ALONG T Hg 5L*NT PATH? 


TOTAL GM7CM»®2 OF WATER VAPOR ALONG THE SLANT PATH, 


IT- IS 


ANStZaj * TOTAL PATH LENGTH JN.CM 


O 9 


00 J I 2 

35® 

THETaS®ThETAS*(-CUN) 

oq» 13 

36® 

PHIL®PHIL*COn 

0011*4 

3 7® 

thetal-thetal* : -con j 

00115 

38® 

DELT 3ABSUE-ZS) 

QO 116 

39 * 

DO 80 I s ! j 32000 

00121 

HU® 

DEkLiDfcLTViD.O 

00 122 

a i ® 

L* I 

00123 

H 2 ® 

IF (DELI «LE .2.0 ) GO ^ Q 

00 125 

H 3 ® 

80 CONTINUE 

QO i 27 

H © 

call exit 



00130 "lb* 6 1 IF (L.LE. 1 ) 0 1 L T »U EL j / 1 0 . 0 

-H01J-2 jttflJt AUS-U LM -4 -t-Q - - ' 

00133 A7* Q1®KE*ZS 

QQ13A *(B« «2*C05lPHlS) 

00135 A9* XS»«1*C0S(THE.TA5)*(k2 

QQ136 SO® . Jt 5.9 W 1 * 5 l N ( T M E T A S ) • (i| 2 „ — . — 

00137 51* HS»Rl*SINtPHIS) 

Q Q1AQ S,2.ft: : - fl 2 ■ CHS If HiLl - 

Q01A1 S3® «1=>RE + ZL 

QQ1A2- .. 5A* AO * QltXQSlTHfcTALl * R 2 . .. - 

00 1 A3 55* YL*Q 1 *S I N ( THET AL ) *<j2 

QQIAA- 56® _ -HO#«l*SlNlPbiU _ _ 

00 l 85 57* A8D“C0SlNV(<lXS*XL) + (Y5*YL)'MHS*HL)>/(5GRT(XS**2*Ys**2«>HS**2) 

00185 5 8* 1»SuRT(XL* »2 +YL* *2+HL *»2 > > ) ' 

00186 59* DO 3 1*1,3 

00151 - 60* .3 Cl I 1*0*0 • 

00151 61* C FROM HERE TO STATEMENT H FINDS THE VECTOR ( C } FROM The TARGET TO THE 

oo.isi 62 * c ... Sa tellite . ....... 1 - 

00153 63* A( 1 ,1 )*SIN{PHIU*C0S(THETAL> 

00 1 5 A 68* A(2,l)»->S1N(THETAL r 

00155 65* A(3»n*C0SlPHlL)*C0SiTHETALl 

00156 66* At 1 j 2 )* SIN (PHlLj • S IJN ( T HE T aL > __ 

00157 67* A ( 2 » 2 1 * C 0 S ( T H E T A L ) 

00160. 66* A( 3i2 1*C0S(PH1L>»S 1 N( THE TA U ... 

00161 69* A ( 1*3 )*~C0S(PHIL) " 

00162 70* A 12, 31*0.0 • 

00163 71* A 1 3 ,3 J *S IN 1PH 1 L ) 

QQ 1 6 A 72* 8(1 ) *XS“XL . _ 

00165 73* 8(2 )*YS“Yl 

00166 23* Si 3 >*J4S"-HL - - -- -• — 

00167 7b* DO A 1*1 ,3 

00172 76* Og A M* 1 , 3 . 

00175 7 7* A Cl J )*A( I »M)*8(M)*C(l 1 

Q0200 78* Ph JL-PHIL/CON __ _ 

00201 79* THETaL*THETAL/ 1-CON 1 

00202 80* PHIS-PHIS/CON 

00203 81* THETaS*TH~£TaS/(-CUN> “ 

0020 A 62* PHI*ATAN2(S0RT (C( l ) **2*C12I**2) ,C(3i 1 

00205 83* IF (PH I ,GT . .01 7) PHl*PHl-«0092833 

00207 8 A » IF 1PHI/C0N.GT .90.0 IftRlTE (6,68) _ 

00212 85* 88 FORMAT (///« IX , ’WARDING, ZENITH ANGLE OF UNREFRaCTEq PATH EXCEEDS 

00212 86*_ 190,0 DEG * . / , 1 X , * IT Is HIGHLY PROBABL E T HAT T HE MftcRAFT QR SPACE 

00212 87* 2 C R A FT CANNOT SEE THE TARGET*,///) 

00213 sa* a<? call modaim ul+uelt*.s,pp „ahalti ,xlamda) 



0021 H 

0 9$ 


phi int»phi 

0011 S 

. 9 Q« 


21*10- - . .. 

0Q216 

9 I « 


D 1 3 ANS £ 3 > 

002 A 1 

92® 


WATER 1»ANS ( S3 ) 

00220 

9 3® 


XN1»ANS(18) 

0.022.1 

9 9 « 


SUK a O * 0 

00222 

9 b $ 


SUmI*0®0 

- 00220 

9 6® — 



-SUM 2^0 « 0 

QQ22H 

9 / ® 


SUM3 s 0<0 


9d« 


Sum9=D.Q . ... _ - - 

00226 

99® 


DO S 1*1 ,32000 

0023 { 

100® 


Z2*n®DE0T. _ 

00232 

SOS* 


SjsREtZl 

Q0233 

102® 


S2«RE®Z2 

GG23R 

103* 


call moOatm (Z2*oeli®.s # pf , hhaot i .xlamda » 

00235 

109® 


02*ANS(3> 

00236 

SOS* 


WATER2>ANS( S3) 

00237 

106* 


XN2 b ANS<18» 

Q02A0 

107® . 


PS I a S s N I NV ( S 1 ®S 1 N S ph I ) /S2 > 

00291 

108® 


PHlPR*55N|NV(Sl*SSNipHl)*XNl/(S2®XN2)J 

QQ292 

l 0 9 • 


0UM B Dl®Sl*SiN<PHI->P5l>/SIN(PS|)®1.0E*a5 

00293 

no* 


SUMi«SUMS>PHI-P5l 

00299 

HI® 


S0M2»SUM2®*ATERS*0UM/l 000.0 

00295 

H 2® 


SUM3®SUM3®0UM 

002*86 

Hi® 


SUM9»SUM9®QUM/U s 

002*87 

i S9# 


IF <22.SE.ZSS_ SO TO 82 

00251 

S i 5® 


SUM“SUM®ABS<PHIPR W PSS ) 

00252 

I S 6® 


phi »ph s pr 

00253 

S S 7 * 


Z | » L Z 

0025*: 

118* 


0 1 ®D2 

00255 

11?* 


WATER 1 ®WATER2 

00256- 

1 2.0® 

i 

XN I ®XN2 

00260 

!2S« 


CALL. EXIT 

00 2 6 L 

■ 112® 

8jL 

continue .. ... •- 

00262 

123* 


Q = 5UM 1 “ABO 

00263 

12*4* 


PHI“PHi iNT-0/2.0 

00269 

125® 


IF ( ABS < q ) *SE. .0001 \ SO TO 89 

00266 

1 26® 


A WS (21 »=*P HI 

00267 

12?* 


ANS <221 «*5UM3 

00270 

128® 


Af+S<23|atSUM2 

00271 

129® 


ANS<29»»SUM9 

r 0 0-2 7 2 

- 4 30* 


if tPHUX0N.Ll.9a. 01 -60 10 -03 - -- - - — -- - - 

002 7 *4 

1 3 S ® 


WRITE <6*878 

00276 

132* 

8? 

format HX.///.SX. * the angle from zenith is greater than ?o.c*j 



ts3> 


00277 

S 33# 

ANS C 22 J “Q « 0 

. Q Q JQ Q 

J 2H • 

ANS 5-23 i»GL.Q 

00 JQ S 

l 3S«s 

ANS(2‘U*0.U 

_ QQ1Q2 _ 

1 3-6 » 

S3 RETURN 

□ 0303 

137* 

End 


ENO OF "UNIVAc lias FORTRAN V COMPILATION* 0 *OIAGNOSTic* MESsAGEOj 



OOSOS 1® FUNCTION COS I NV 4 A 9 

00 1 0 1 'l9 C THIS FUNCTION calculates the INVERSE cosine of *a*« 
QO { 03 3® CQS1NV«ATAN2 (S«KT ( i »0-A**2 J ,-A > 

QQIQH .. ^ « RETURN .! 

00105 5® END 


END OF UNIVAC UGH FORTRAN V COMPILATION, 0 aOlAGNQSHC® NfcSsAGECSj 


00101 
00101 
00 I 03 
OfllQ^ 

00105 


1 • FUNCTION SININV(A) 

2® C THIS FUNCTI ON_C ALCUL A TE S THE INVERSE SINE OF »A*, 
3* S I N 1 N V <■ A T A N 2 (A, 4 S 0 R l • 0 ** A * ® 2 1 ) 1 

4* RE TUR N • ■ : 

5 * END 


end OF UNIVAC 1108 FORTRAN V COMPILATION, V sDlAGNOSTIC® MESsAGE(S) 


jamjQJ UL FUNCTION J LLP-tli 

00101 2* C W ■ SPECIFIC HUMIDITY to I TH UNITS OF GM/KG 

00 IQ 1 __3 • C SPECIFIC MUMlDlTYaGM OFftATER VAPOR / 4*G OF a I R . I NCL Up I NG to ^TER V APO R) 

00103 *4® X*E(T> 

00104 5® 0*Q»621 < ?7®X/ (P~0«37Sq3®X J « J 0 0 0 • 0 1 __ 

00105 6* RETURN 

00106 7® END 


END OF UNIVAC 1106 FORTRAN V COMPILATION. 0 ^DIAGNOST IC® MESSAGE ( S ) 





00101 

i * 


00 1 0 1 

2 * . 

.fL.. 

00101 

3 » 

Q * ® W « 

0Q201 


C 

00201 

b* 

C G 1 V 

00102 

6# 



001 Li 
00112 
00111 . 
OO H <4 
00115 


0012 6 


FUNCTION ALTJTU lTVMI(,h»TVUOiH,PHiQH*PLO«,HUOfll 


I ® * -8 


:n the: temperature and pressure at each of z points and The altitude qf 
the lower point, t h is fu noti on calculates the altitude of the higher point 

A CONSTANT « -m*g/r 


00 1 0 I 

7* 

C 

ALTITU IS IN METERS* CONN IS 

QQ1Q.1 

_ 8* 

C 


00IC1 

9® 


□ 0102 

10* 

C 


00 103 

l 1* 


DATA C0NN/-3.R16319‘47E”02/ 

00105 

12# 


D*T VH IGH-TVLOW 

00106 

1 3 * 


1F(D) 2.3*2 


I • C I 


IN* 
l 5® 
1 6 » 


17* 

18 m 


2 D«_ ; _ ■ 

ALTITU ■HLOW+ITVHI&R-.TVLOWI/D 

GO TO 6 _ 

3 ALTITU »Hi.UW*TVLOW*ALOG(PHiGH/PLOW)/CONN 

6 RETURN . ' 


19 # 


end 


END OF UNIVAC 1108 FORTRAN 9 COMPILATION. 


o .diagnostic* message is » 



OOIOI 1® FUNCTION >REs'lPLo¥ 8 ^'TVLoi®TVHlSH»OHI ‘ 

00103 2* DATA C0NN/-3-. <4 1 63 1 i*<7E.-Q2/ . 

00103 3* C . 

OQ 1 03 ... 

00103 5® C . ' . ' ' ' ~ 

00103 6 ® C THIS PROGRAM CALCULATES PRESSURE -PRES- AT SOME POINT -PH- ABOVE A 

00103 7® C POINT IN THE ATMQSPH E RE HAVING PRESSURE -PLO*- »HEr£ -Q- IS THE 

00103... 8 ®, .„_£_' ; TEMPE RATU RE <*N A 0 J E N j. ft NO. -TVHIGH- and -TVL OW" A«E CORRESPON D IN G 

00103 V® C TEMPERATURES. -conn- IS CONSTANT ■ -M*G/R 

Q. 010.3 t o® c ‘ ' _ _■ ; 

00 103 II® .C**®®*®®*®®®®®*®*.®***®**®®®®*®**®*®'®*®*®*®®®®*®®®***®*®®**®*®*®®*®®*®®**®®®®*®** 

00103 12® C ^ : 

00105 S3® IF ( D ) 2.3.2 

001 10 l_H* 2 PRES*PEO** { TVHI 6H/TVU0A ) ** (CONN/D ) ' . . 

001 11"“-"“ 15* V GO TO *4 ' ~~ r ~ " ~ ' ‘ 

00.11.2 .._|A® 3 PKES»P t.O» *£XP <-CQNN*PH/TVLOW ) __ _ _ 

00113 17® *4 RETURN ~~~ ~ ' ~~ : ~ 

001H 18 ® END 


END OF UN I VAC 1108 FORTRAN V COMPILATION* 0 «0 1 AGNOST l C* HES gAGElSl 


cn 


a 0.1.01. 

00103 

00 103 

00103 

00103 


. 1 * 

2 * 

3 ± 

*4® 

5* 


FUNC TIO N £ t X I 

DATA TS/373. 16/.T0/273.16/ 


L 




c 


00103 6® 

00 LQ 3 7® 

00103 8® 

0 0103 9® 

00103 10® 

001 03 11® 


c This ROUTINE calculates VAPOR pressure OVER a plane surface of 
c I A TER CC„*._0« L g» OR...UF. ICE.. 273.16 ) BASED ON TEmPERaTMFE IN DEG 

C KELVIN. E i X ) IS IN MB ' 

C._ SEX. C ® 2 7 3 » . 1 6 IF . Y .0 U u A NT__ V A P 0 RjP RES. OV ER ICE USED BELO m, 273. DE G K 

C ' 


»«••••< 


00103 12® C 

QflLOA i.3% Q ._ _ 

00107 1*4® T = X-C 

HOLLO 15-t _ l£__lX-®l.£®._Jl*Q-l. -.60- _L.Q_.*1 _.... - 

00112 18* IF ( T ) 1 ,2,2 

00112 17® C FORMULA FOR VAPOR PRESSURE OVER IC E 

00115 | 8* I E»6,107i*10.0®*l-9.09 718* ( - l . 0 + T0/ X J -3 . 566S‘4®L06 l 0 I TO/X ) *CK P W ' 

'_Q.Hl 1.5.-_ _L.il 1. .• U--«..0.“ X./XO..I .L .. ". . V ' : 

00116 20 * 60 10 5 

00 J 1 6 2 C _ . _ __ 

00116 2 2 9 £*©©*©*•*©©©©•©•©*♦•©©♦*•**♦©•©♦*♦•♦*©*©©*©*©♦♦©♦ ©©©©ifc©©©®©*©©©®©©©©©*®®©#**^^,^®^ 

QC1I6 23® C 




os 


0011 6 2H9 C FORMULA FOR VAPOR PRESSURE “OVER WATER 

00117 2b* 2 £-i0i3 e 2 t )6»i0 s Q**i-7i, 90 298»t-»l t 0-t-TS^Ai-»5tU 2aue*L06lUITs/X>"i«38l6E 

001 17 20* i-U7®UG*G SMI, tn«'3‘lH«>1|,0-x/TSn“l®Ql*a<>l328E“U3®U0,0® l M“3® 4 49ji‘«*!'"l 

00147 27* 2®0 *TS/X t 1«1»0) 1 • ; 1_ 1 

00120 28* GO TO 5 

0012 1 29* Jf £ = 0 « G ... _ ; ■ . ... ' 

00122 30* 5 return 

00123 31* ENQ ... . _ i_ 

END OF UNIVa C 1 108 FORTRAN V. COMPi (.A T I 0 Nj 0_ » 0 1 aJLN OSTlc. H£SsaG£ ( S; 


001 04 - 

00101 
0Q101 
GO IP 1 
00 104 
00101 
00101 

00 1 0 i 

001 01. 
00101 
QO IQS 
00101 
00 1 03 

00106 

00107 

00107 

00 1 l 0 

00111 

00112 

001 13 


JjL 


J-.imCX 40 LN._ RJL SU E_*JL1 . 


2 * 
3 * 
<1* 
4ml. 

6 * 

7 » 


x this routine calculates the mining ratio (gm of.-h^oi/uo of pRy 
c based on x which is temperature in oeg kelvin 
c KlSiP.m *0/00 (IE PARTS PER THOUSAND) 


8* 

9 * 

10* 

4L* 


12* 

13 * 


C S IS VAPOR PRESSURE OP WATER 

c p _is total atmospheric pressur e in Mg 

C 


IF IS) 7,6,7 


1*1® 7 CONTINUE 

I S* _ R*13.Q16»S*F(P,X)/U8«966 < t»tp-S«F(P»X) ) 1* 1 000. 0 

16 * C R IS IN GM/KG 

l 7 » RETURN __1 

1 B* 6 R*O.Q 

„li* RETURN _____ 


20 * 


END 


END OF UNIVAC 11Q8 FORTRAN V COMPILATION, 


0 *OlAGNOSTIC* MESSAGE f S 1 



00101 
00 103 
0 0 1 0 *1 
0 0 I 0 s 
00 i OS 

00 1 os 

GOiOS 

00 1 os 

QQSOS 
00 1 os 
00 i OS 

00 10 s 
00 1 os 

00 l OS 
0010S 

00 1 OS 
00 1 OS 
OQ 1 QS 
00 1 OS 
00 1 OS 
00 1 OS 

oo no 
oo n i 
oo n s 
00116 
00 117 
00121 
00 122 
00 123 
00126 

00130 

00131 
00133 

00 1 3S 
0 0 13 5 
00 136 
QO 1 37 
00 1 SO 

o a i s i 

0 0 1 S 2 
00 I S3 
00 I SS 
0 0 1 S 5 


1 * 

2 

j » 
H * 

b» 
6 $ 
7 « 
& * 
9 * 
1 0® 
1 1 * 
t 2 ® 

1 3 » 
1 S * 
lb* 
1 6® 
17® 
16® 

1 V a 
20 ® 

2 1 » 
22* 
23® 
2 S « 
2b* 
26® 
2/« 
2b® 

2 7* 

3 U * 
3 i * 
32 * 
3 3* 

jl* 
3 b * 
3 6 * 
3 7 * 
3 8 * 
3 9® 
HU® 
HI® 
H 2 • 
H 3 ® 


function f ( p f x j 

DIMENSION T£( i2*»FE U i l»U( U 3 U ) 

U A T « llU(lfJ) »J = lsi i )si ss i8|2^ /Q#»i* jS3»s 

226© ® 3 6 ® g H 6 ® j bb ® 5 i© §2© s 3 ® s H © s 6 o 9 I 1 « $18® &2H®&3H» e 3^ ® ? I « 8 2 ® iH© © 

38® ^7® g I 1 ® s i b © ? ® §32® t H l 9 $ H9 e> * 0 © & 2 ® s 5 ® 56 ® $ 8 © ? J, 2 ® §16®®2H© $32® gHO® & 

H H 7 ® $ H $ Q ® 2 1 U® 9 IS® 9 i 3 ® j © g 3 2 9 § H (J « 9 M 7 ® g H ** 0 ® 9 12® 9 1 & ® ® 2^ ® * Z~? © & 3H® rl i ® 1 
b H 8 © s23« s JO ® j i / 8 s HI e 5 5 Q ® §6^0© g 2 3 © g 3 H ® 9 H l ® g H B ® & S ^ ©7^0® 8 3 / © **$£>• 

O * b 2 ® s b 9 ® » 8 ^ Q © » ^8 © sbo® ® 6 H ® / 3 f ^ 5 Q ® j *40® 9 3 U * 0 *® 2 0 ® 0 « i L) ® g 0© 9 10® & 2 Q ® 9 

7 3 U • ^40° e b U ® 3 6 0 ® 7 9 P E ^ S © 0 1 G © g i(J o & 5 0 ® ? 1 G 0 ® g 2 Q U © g 3 U U © s bOU © g 7 Q G ® & 9 Q Q @ & 

B i 1 UO ® 7 • _ 

C 

C 

C *F' 1b The CORRECTION FACTOR FOR THE L>fcH M KTUR£ OF T ME m 1 A I u K £ OF MR 
c a ij [j water vapor from the i u f. a l gas i a « # 

C A IS TEMPERATURE in DEG K t l v. I n 
0 P lb Total A T MO b PHE R 1 C PRESSURE IN mb 
c 

C*»®®**^#«A#**»»»*»**»aa*»*»*****»«***<ft*®#»e*****S»®***S®»**»*i0««®®««# 

c 

T=X-273»16 
00 l 1=1,12 

IF 1 T . LE * T E i I ) ) GU f o 2 
11*1 

1 CONTINUE 
F A = * .0 

0 ' f 0 ^ 

2 00 S J= 1 , 1 1 
IF (PsLtaPESJ) J GO ^ 0 b 
J U * U 

H CONTINUE 
F A * 1 » U 
GO TO 3 
b I*ll 
J = JJ 

V 1 * 1 U 1 l*! i J i ”U ( 1 s J ) )/l0»0®{T-TE( 1 t l + 0( 1 « J 1 
F 2 =(Uli + l 80 +l)-Oll»b+i))/io» 0 »(T-TEti))+UC 1 . 0 +li 
FAalF2«Fl)/(Ft(J+l)Ht(J)MlF-PElJI»+Fl 
F A = 1 • U + F A * 1 eOE-OS 

3 F = F A 
RETURN 

end 


E_NO. OF UNlVAC UQ3 FORTRAN v COMPILATION 


0 * 0 J A G N 0 S T I C * H £ 5 a a g E 1_S 1 



